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Arcsecond-Resolution Submillimeter HCN Imaging of 
the Binary Protostar IRAS 16293-2422 
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ABSTRACT 

With the Submillimeter Array (SMA) we have made high angular-resolution 
(~ 1"= 160 AU) observations of the protobinary system IRAS 16293-2422 in the 
J = 4 - 3 lines of HCN and HC^^N, and in the continuum at 354.5 GHz. The 
HCN (4-3) line was also observed using the JCMT to supply missing short spacing 
information. Submillimeter continuum emission is detected from the individual 
binary components of Source A in the South-East and Source B in the North- West 
with a separation of ~5". The optically-thin HC^^N (4-3) emission taken with 
the SMA has revealed a compact (~ 500 AU) flattened structure (P. A. = -16°) 
associated with Source A. This compact structure shows a velocity gradient along 
the projected minor axis, which can be interpreted as an infalling gas motion. 
Our HCN image including the short-spacing information shows an extended (~ 
3000 AU) circumbinary envelope as well as the compact structure associated 
with Source A, although the P.A. of the compact structure (= -45°) seen in 
the HCN emission is slightly different from that of the HC^^N emission. A toy 
model consisting of a flattened structure with radial infall towards a 1 Mq central 
star reproduces the HCN/HC^^N position- velocity diagram along the minor axis 
of the HC^^N emission. In the extended envelope there is also a North-East 
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(Blue) to South- West (Red) velocity gradient across the binary ahgnment, which 
is hkely to reflect gas motion in the swept-up dense gas associated with the 
molecular outflow from Source A. At Source B, there is only a weak compact 
structure with much narrower line widths (~ 2 kms~^) seen in the optically- 
thin HC^^N emission than that at Source A (> 10 kms~^), and there is no 
clearly defined bipolar molecular outfiow associated with Source B. These results 
imply the different evolutionary stages between Source A and B in the common 
circumbinary envelope. Our study demonstrates the importance of adding short 
spacing data to interferometer data in order to probe the detailed structure and 
kinematics from extended (> 3000 AU) envelopes to inner compact (< 500 AU) 
structures around low-mass protostars. 

Subject headings: ISM: individual (IRAS 16293-2422) — ISM: molecules — stars: 
disk — stars: formation 



1. Introduction 



Millimeter molecular-line observations have found 2000-10000 AU scale molecular en- 



velopes around low-mass protostars 


Benson k Mvers 1989: 


Zhou et al. 19891: 


Blake et al. 1994 




Zhou et al. 1994: 


Blake et al. 1995: 


Hogerheijde et al. 1997: 


Hoeerheijde et al. 1999: 


Takakuwa et al. 2000 



rotating and infalliii g gas motion in these envelopes (lOhashi et al. 19961 : lOhashi et al. 1997a 



Ohashi et al. 1997bl : iMomose et al. 19981 ). The dense and warm innermost (< 500 AU) 



region of these low-mass proto stellar envelopes is a likely site of formation o f protoplan 



etary disks around protostars (ITerebey et al. 19931 : iBeckwith fc Sargent 19961 ). However 



there is not a great deal of information published about the connection between large-scale 
protostellar envelopes and small-scale disks around the central protostars. With the ad- 
vent of interferometers in the millimeter and centimeter wavelengths and the use of mul- 
tiple spectral lines, progress is being made toward separating the warm and dense regions 



from the overlying lower-density and colder gas along the line of sight (IMundv et al. 1990 



Momose et al. 19981 : IPuUer fc Wootten 200ol : iHogerheiide 200ll : iTakahashi et al. 2006h . Re- 
cent arcsecond or subarcsecond-resolution millimeter interferometric o bservations have en- 
abled us to investigate the innermost reg ions of protostellar envelopes (ILooney et al. 2000l : 



Beltran et al. 2004J : iBottinelli et al. 20041 ). Detailed submillimeter line observations are also 
useful for these studies, since submillimeter molecular lines such as HCN (4-3) trace higher 
densities (> 10^ cm~^) and temperatures (> 40 K). Earlier studies of low-mass protostel- 
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lar envelopes with the Submilhrneter Array (SMA)Q have demonstrated that submilhrne- 
ter interferometric observati ons are uniquely suited to probe the i n nermost part of low- 
mass protostellar envelopes (Kuan et al. 20M : Ixakakuwa et al. 2004 : Chandler et al. 2005 : 
Bourke et al. 20051 : |j0rgensen et al. 20051 ). We note that in order to study the connection 



between the large-scale envelopes and inner disks, observations that sam ple adequately the 
emission from large to small scales with h i gh enough resolution a re needed ( 
Welch et al. 200ol : IXakakuwa et al. 20041 : IXakahashi et al. 2006h . 



Gueth et al. 1996 



IRAS 16293-2422 (hereafter 116293) is a Class protobiiiary system with a projected sep- 

aration of ~ 800 AU (Source A and Source B) JWootten 19891 : iMundv et al. 19921 : lloonev et al. 200oh 
i n the Ophiuchus Molecular Cloud, sur rounded by an ~ 8000 AU-scale circumbinary envelope 
(ISchoier et al. 20021 : IStark et al. 20041 ). Previous SMA observations reported by Chandler et 
al. (2005) further demonstrate that one of the binary components. Source A, may be com- 
posed of at least three components within a 1" region. It has been re ported that there is 



rotating and infalling gas motion in the circumbinary envelope of 116293 (IWal 



Menten et al. 19871 : IZhou 19951 : iNaravanan et al. 19981 : ICeccarelh et al. 2000al: ISchoier et al 20021) . 



ker et al. 1986 



Large-scale quadrupolar molecular outflows have been obseryed in 116293 (I Walke r et al. 1988 



Mizuno et al. 1990: Castets et al. 2001 



Hirano et al. 2001 



Garav et al. 2002: iStark et al. 2004J). 



Previous interferom etric ( Bottinelli et al. 2004 : Kuan et al. 20041: Schoier et al. 20041: 



Huang et al. 20051) and single-dish molecular-line observations ( 



Ceccarelli et al. 2000b: Schoier et al. 2002: Cazaux et al. 2003 



Blake et al. 1994 



Chandler et al. 2005 



van Dishoeck et al. 1995 



Schoier et al. 20041 ) have re 



vealed a number of complex organic molecules and hot-core chemistry in 116293. Previous 
published SMA observations (IChandler et al. 20051 ) have studied in detail the structure and 
kinematics in the circumbinary envelope and circumstellar disks for each binary companion, 
the driving sources of the outflows, and the evolutionary status of each member of the binary. 

In this paper, we describe results of HCN (4-3), HC^^N (4-3), and 354.5 GHz continuum 
observations of 116293 made with the SMA and James Clerk Maxwell Telescope (JCMT). 
This paper focuses on sampling, as completely as we can, the structure and kinematics 
from the large-scale (~ 3000 AU) circumbinary envelope traced by JCMT to the smaller- 
scale (< 500 AU) structures traced with the SMA. For comparison wi th previous studies 
we assume a distance of 160 pc to the Ophiuchus cloud (IWhittet 1974r). but we note that 



more recent studies locate the cloud at a nearer distance of 120-140 pc (Ide Geus et al. 1989 
Knude fc H0g 1998h . 



^The Submillimeter Array (SMA) is a joint project between the Smithsonian Astrophysical Observa- 
tory and the Academia Sinica Institute of Astronomy and Astrophysics, and is funded by the Smithsonian 
Institution and the Academia Sinica. 
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2. Observations 

2.1. Submillimeter Array Observations 

With the Submilhmeter Array (SMA) we made HCN (J = 4-3; 354.5055 GHz), HC^^N 
(J = 4-3; 344.200122 GHz), and 354.5 GHz continuum observations of 116293 on 2003 
March 14, July 12, and 2004 June 19. Details of the SMA are described by Ho, Moran, & 
Lo (2004). The SMA is a double-sideband instrument, and the HCN and HC^^N lines were 
observed simultaneously in different sidebands. We assigned a spectral window of the SMA 
correlator ("chunk") with a spectral resolution of 0.17 kms~^ to the HCN line, and a chunk 
with a resolution of 0.71 kms~^ to the HC^^N lin^. Table 1 summarizes the observational 
parameters. The observations were made in three different array configurations to provide 
well-sampled {u, v) coverage. The range of the baseline length projected on the sky was from 
~ 10 kA to ~ 223 kA, and our observation s were insensitive to s tructures more extended 



than ~ 16'.'5 (~ 2600 AU) at the 10% level flWilner fc Welch 1994[ ). We confirmed that the 



visibility amplitudes of the continuum emission from 116293 from the three observing periods 
were consistent within the noise level. The overall flux uncertainty was estimated to be ~ 30 
%. Part of the 2003 data set has already been published by Kuan et al. (2004) and Huang 
et al. (2005), which is of lower angular resolution than the full data set presented here. 

The raw visibility data were calibrated and flagged with MIR, which is an IDL-based re- 
duction package adopted for the SMA fro m the MMA softwa re package developed originally 



for the Owens Valley Radio Observatory (jScoville et al. 19931 ) . The calibr ated visibility data 



were Fourier-transformed and CLEANed with MIRIAD to produce images (jSault et al. 19951 ), 
The spatial resolution is I'.'l x 0'.'6 (RA. = 39°), 1'.'3 x 1'.'2 (RA. = 30°), and 1'.'6 x 1'.'3 (RA. 
= 14°), in the 354.5 GHz continuum, HCN, and HC^^N images, respectively. The spatial 
resolution of the HCN image is set to be the same as that of the combined SMA + JCMT 
image (see next section), in order to directly compare the SMA and combined images and 
to see the effect of missing short-spacings. 



2.2. James Clerk Maxwell Telescope Observations 

Submillimeter single-dish mapping observations of 116293 in the HCN (4-3) line were 
made with the James Clerk Maxwell Telescope (JCMT) on 2004 March 7 and April 7. The 
JCMT beam size at 354 GHz band was ~ 15" and the typical system temperature was 600 



^Due to the instrumental constraint of the SMA correlator, it is not possible to asign the high-resolution 
chunk to both the HCN and HC^^N lines. 
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K. We observed a 7 x 7 map centered on the field center of tlie SMA observations at a grid 
spacing of 7'.'5, providing a Nyquist-sampled map in the 45" x 45" region. Each position in 
the map was observed for 1 minute and the map was repeated 6 times. The central map 
position was observed at the start, middle and end of each map as a check on the relative flux 
calibration. The observations were made in dual-polarization mode, and the 2 polarizations 
averaged, resulting in an rms noise level per 0.13 km s~^ channel of ~ 0.13 K in T^. Pointing 
and focusing were checked before each set of maps. The conversion factor from T^ (K) to S 
(Jy beam~^) was derived to be 36.7 as; 

^ Vmb 

where ks is the Boltzmann constant, ^Ibeam is the solid angle of the JCMT beam (= 
15"), A is the wavelength, and rjmb is the main beam efficiency of the JCMT (= 0.63). The 
mapped region covers well the field of view of the SMA observations (~ 35"), which enables 
us to compare the JCMT fiux to the SMA fiux. The SMA observations recovered ~ 35% 
(Source A) and 31% (Source B) of the total HCN(4-3) fiux observed with the JCMT toward 
the binary system (The JCMT observations do not resolve the 5" binary). We combined 
the JCMT data with the SMA data, and in the subsequent sections we will mainly discuss 
the combined images. The details of the combining process are described in Appendix. The 
resultant synthesized beam size and the rms noise level in the combined images is 1"3 x 1"2 
(P. A. = 30°) and ~ 0.90 Jy beam~^ channel"^, respectively, where the velocity resolution is 
the same as that of the SMA data. 



3. Results 

3.1. Submillimeter Continuum Emission 

In Figure 1 we show our 354.5 GHz continuum image of 116293. Previous SMA im- 
ages of the continuum emission in 116293 at somewhat lower angular resolutions have been 
presented in Kuan et al. (2004) and Chandler et al. (2005; although they subsequently 
improved their resolution through "super- resolution" imaging techniques, discussed below). 
Submillimeter continuum emissio n from the individu al binary protostars, called Source A 



(south-east) and B (north-west) (IMundy et al. 19921 ). are evident. The total continuum 



fiuxes at 354.5 GHz are 3840 and 4050 mJy at Source A and B, respectively, which are 
slightly higher than the continuum fiuxes at 305 GHz (3460 and 3150 mJy at Source A and 
B; Chandler et al. 2005). Chandler et al. (2005) reported that the spectral indexes (= a) 
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in the millimeter and submillimeter regime are 2.91 and 2.51 at Source A and B, respec- 
tively. We re-estimated the spectral indexes including our new submillimeter measurements, 
and the values are 2.87 and 2.68 at Source A and Source B, respectively. These values are 
consistent with the values by Chandler et al. (2005) within the uncertainties (~ ± 0.1). 
The submillimeter continuum emission at Source A shows evidence of extension along the 
North-East to South- West direction, while that at Source B is more compact. In Table 2, 
we summarize the properties of the submillimeter continuum emission from both protostars 
after the deconvolution of the synthesized beam. Using "super-resolution" imaging by only 
using data with uv distance greater than 55 fcA, the extension at Source A is resolved into 
two components with the SMA (Aa and Ab; Chandler et al. 2005), possibly a close binary 
system. The peak posi tions of the 354.5 G Hz continuum emission are slightly (~ 0'.'4) off- 
set from the 115 GHz Jloonev et al. 200oh and 300 GHz Jchandler et al. 2005h continuum 
positions, presumably due to the limited calibration accuracy of our higher-frequency obser- 
vations. Detailed multi-wavelength (230-690 GHz) analysis of the submillimeter continuum 
emission will be presented in a forthcoming paper. 



3.2. Spatial and Velocity Distribution of the HCN and HC^^N Emission 

In Figure 2, we compare the HCN (4-3) line profiles from SMA and combined SMA+JCMT 
data at Source A and B, along with the SMA HC^^N (4-3) line profiles. We also show 
the average (Source A + Source B) SMA and SMA+JCMT spectra for comparison with 
the JCMT HCN spectrum. The HCN (4-3) line profiles show two emission peaks with a 
brighter blue-shifted peak and an absorption dip. The SMA and combined HCN spectra 
show a larger brightness temperature with more prominent wing-like emission as compared 
to those obtained with the JCMT, while the interferometric-only spectra systematically miss 
lower- velocity emission around Vlsr = 2-4 km s~^ as compared to the combined spectra, 
particularly at Source A. The dip in the SMA+JCMT and the SMA spectra seems to show 
absorptions against the continuum. The HC^^N spectra shown in Figure 2 include new data 
taken on June 19, 2004 as compared to the HC^^N spectra shown by Kuan et al. (2004), 
and the spatial resolution (~ 1'.'6 x 1'.'3) is higher than that of Kuan et al. (2004) (~ 
2'.'7 X 1'.'3). At Source A, the HC^^N hue profile shows a single peak near the dip veloc- 
ity. The centroid velocity of the HC^^N line at Source A is estimated to be ~ 2.9 km s~^ 
from a single-component Gaussian fitting to the spectrum. At Source B, there is a weak 
(~ 4cr) HC^^N emission near the centroid velocity. This velocity is slightly bluer than the 
HCN dip velocity (~ 4.3 km s~^), although the coarser velocity resolution of the HC^^N 
data (~ 0.71 km s~^) prevents us from making the direct comparison between the HCN 
and HC^^N velocities. Taking into account the coarse velocity resolution in the HC^^N 
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line, we adopt the symmetric velocity of the HCN Position- Velocity diagram (Figure 7) 
as a systemic velocity (~ 3.6 km s^^). At this velocity, the HCN (SMA+JCMT)/HC^^N 
(SMA) line intensity ratio toward Source A is ~ 4, which implies an optical depth o f 
~ 0.3 for the HC^^N line assuming an N^^/N^^ isotopic ratio of 270 Jhucas fc Liszt 19981 ). 
Therefore, the HC^^N line is likely to be optically thin. The redshifted self-absorption in 
the HCN spectra against the contin uum at So u rce A could be a sign of the inward mo- 



tion toward the continuum source (Zhou 



Di Francesco et al. 2001: Schoier et al. 2002 



992 



Myers et al. 19951 : [Narayanan et al. 1998 



Stark et al. 20041 ) , although the presence of the 



extended HCN absorption in the entire region (See Figure 4) implies significant contami- 
nation from extended low-temperature material. The red-shifted self-absorption from the 
syste mic velocity can be reprodu ced using models of infalling cores with the Larson- Penston 
flow (ILarson 19691 : iPenston 19691 ). as has been shown by Masunaga fc Inutsuka (20 00). Such 
a red-shifted dip was also detected in B335 (IChoi et al. 19991 : lEvans et al. 20051 ). which is 
one of the most well-studied infalling cores. 

Figure 3 shows two different HCN integrated intensity maps; one made using only the 
SMA data (left), and the other made using the JCMT and SMA combined data (middle). 
Figure 3 also shows the HC^^N integrated intensity map taken with the SMA for comparison 
(right). The continuum emission has been subtracted before forming these maps. In the 
optically-thin HC^^N emission there is a compact emission associated with Source A. From 
a 2-dimensional Gaussian fitting to the image, the FWHM size of this compact HC^^N 
emission after the deconvolution of the synthesized beam is 450 AU x 250 AU (P.A. = 
-16°). A similar compact structure is also seen in the SMA and combined SMA+JCMT 
HCN images with a size of 630 AU x 320 AU (P.A. = -43°) and 740 AU x 440 AU (P.A. 
= -45°), respectively, although the P.A. of the HCN component is slightly different from 
that of the HC^^N emission. The HCN emission also shows extensions to the North- West 
and to the North-East of Source A in addition to this compact structure. Importantly, the 
combined HCN image shows that all of these features are surrounded by a halo component 
with a size of ~ 3000 AU, demonstrating that the SMA is not sensitive to the extended 
structure. This halo structure is most likely a molecular envelope surrounding the 116293 
protobinary system. Similar extended structures are also seen in the H2CO (4i 3-31^2) and 
SO (77 '-6b) emission, althou gh the detailed distributions are somewhat different from each 
other Jchandler et al. 20051 ). 



In Figure 4, we present velocity channel maps of the SMA-I-JCMT HCN data. At high 
blueshifted velocities (—2.8 kms^^ < Vlsr < -1-2 kms~^), there is a compact gas compo- 
nent associated with Source A. This high-velocity component corresponds to the compact 
structure seen in the total integrated intensity map of Figure 3. From Vlsr = -0.73 kms~^, 
a secondary component is seen at the south-east of Source B. From Vlsr ~ 1-9 kms~^, an 
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extended halo component appears until Vlsr = 4:A kms~^, where little emission is evident 
due to the extended self-absorption in the HCN emission. After this velocity the redshifted 
halo component appears again. From Vlsr ~ 6.5 kms~^a compact gas component toward 
Source A is seen, which corresponds to the redshifted counterpart of the compact structure 
at Source A. 

In order for us to see the systematic velocity structures more clearly, we integrated these 
velocity channel maps into four different velocity regimes, that is, low-velocity blueshifted 
(2.5 - 3.4 km s~^) and redshifted (5.0 - 6.5 km s~^), and high-velocity blueshifted (-2.6 - 
0.7 km s~^) and redshifted (7.4 - 9.1 km s~^) emission. In Figure 5, we co mpare these fou r 



HCN velocity channel maps to the CO (2-1) maps taken with the SMA ( lYeh et al. 20071 ). 
integrated into the same four velocity regimes. The spatial resolution in the CO (2-1) map 
is 3'.'3 X 2'.'0 (P. A. = 40°). In the low- velocity regimes, the overall distributions of the HCN 
emission, which shows extended halo structures, are very similar to those of the CO (2- 
1) emission. This suggests that the extended halo structures are affected by a low-velocity 
outflow or turbulent gas traced by CO. The redshifted halo component appears mostly at the 
south-west and south-east of the binary, while the blueshifted halo component at the north- 
east and north-west of the binary. In the high velocity regimes, the HCN emission shows 
extensions as well as the compact structures associated with Source A. The comparison 
with the CO maps shows that the HCN extensions have distributions similar to those of 
CO, particularly at the high redshifted velocity. This suggests that the HCN extensions are 
also affected by the high-velocity outflow. On the other hand, the compact HCN emission 
associated with Source A seems to show different distributions from those of CO, which 
suggests that it has different origin. The compact HCN component is also distinct from the 
more extended low-velocity emission. We will discuss the origin of these different velocity 
components in the next section. 



4. Discussion 

4.1. Origin of the Different Velocity Structures 

Our SMA and JCMT observations have revealed a high-velocity compact (~ 500 AU) 
structure associated with Source A and an extended (~ 3000 AU) low-velocity circumbinary 
envelope around sources A and B. In this subsection, the origin of these different velocity 
structures is discussed. 

In Figure 6, we compare the high-velocity blueshifted and redshifted HCN and HC^^N 
emission to the CO (2-1) outflow map taken with the SMA (Yeh et al. 2007). As discussed 
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in Yeh et al. (2007), source A is the most likely to be driving the E-W CO outflow, while 
the driving source of the NW-SE outflow remains uncertain. In the optically-thin HC^^N 
emission, the high-velocity blueshifted and redshifted components exhibit a velocity gradient 
along the minor axis (P.A. 74°) of the compact HC^^N structure. The velocity gradient along 
the minor axis is approximately parallel to the E-W CO outflow, though the position angle 
of the E-W CO outflow can not be well-deflned due to the wide opening angle. Although the 
direction of the minor axis of the compact HCN emission is slightly different from that of 
the compact HC^^N emission, the direction of the velocity gradient in the HCN emission is 
similar to that of the HC^^N emission. One possible interpretation of the observed velocity 
gradient in the compact flattened structure is outflowing gas, because both the CO outflow 
and the HCN/HC^^N emission show the same trend of velocity gradient. In fact, the extended 
HCN emission is strongly affected by the molecular outflow, as has been discussed in the 
previous section. The slight difference of the position angle between the compact HCN 
and HC^^N component is also likely to be due to the contamination from the extended 
component. However, the optically-thin HC^^N emission, which should be less affected 
by the extended outflow component, is clearly elongated in a different direction from the 
outflow direction (Figure 5 and 6), suggesting that the velocity gradient is unlikely to be 
due to the outflow. We interpret the observed velocity gradients in the HCN and HC^^N 
emission to be an infalling gas motion i n the flattened diskl ike structure toward Source A 
flHavashi et al. 19931 : lohashi et al. 19961 : iMomose et al. 1998h . 

In the left panel of Figure 7, we show position-velocity (P-V) diagrams along the minor 
axis of the HC^^N compact emission at Source A (see Figure 6), in the HCN emission with 
the SMA-only (gray) and SMA+JCMT (black contour), and in the HC^^N emission with 
the SMA (red contour). In the P-V diagrams, both the HCN and HC^^N lines trace the 
higher-velocity emission located close to Source A, and this velocity structure corresponds 
to the compact emission at Source A. We made a number of model P-V diagrams of a 
geometrically thin disk with Gaussian intensity distribution to reproduce this component, 
one of which is shown in the right panel of Figure 7. From the estimated major and minor 
axis in the SMA HC^^N image we estimated the inclination angle of the flattened structure 
from the plane of the sky to be ~ 57°(= cos~^(maj/min)), and we adopted this inclination 
in the mode l P-V diagram. W e also assumed the internal velocity dispersion (= a) to be 
1.0 km s~^ ( IStark et al. 20041 ). On these assumptions, we made model P-V diagrams with 
different central masses, and the model P-V diagram in Figure 7 shows that the compact 
high-velocity emission could be interpreted as an infalling disk onto the central stellar mass 
of ~ 1 M0 with an acceptable mass range from ~ 0.5 Mq to ~ 2.0 Mq. An uncertainty 
of ±10°of the inclination angle of the flattened structure could reproduce additional ~ 20% 
uncertainty of the central mass. From VLBI studies of H2O masers in 116293, Imai et al. 
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(1999) have found a rotating-infalling disk with an outer radius of 100 AU around Source 
A. Their estimates of the inchnation angle and the central stellar mass are ~ 35°and ~ 0.3 
Mq, respectively. These estimates are roughly consistent with our estimates (~ 57°, and ~ 
1.0 Mq), and the compact HCN structure detected with our observations can be interpreted 
as tracing the outer part of the same infalling disk. 

From the optically-thin HC^^N emission, we estimated the mass of the compact structure 
within the radius of ~ 230 AU (= r) to be 0.084 Mq {= Mr). Here, w e assume the HC^^ 
emission to be optically thin, the excitation ter nperature to be 3 (K) (ISchoier et al. 20021 ). 
and the HC^^N abundance to be 7.4 x 10"^^ (IKuan et al. 20041 ). The infaUing velocity at 
that radius is 2.8 km s~^ (= Vr) from our model. Then, the mass accretion rate (= M) can 
be estimated from 



M = MrVr/r (2) 

and we find a value of 2.2 x 10~^ (M© year~^). This value of the mass accretion rate 
is 4 times higher than that estimated by Schoier et al. (2002) (~ 5 x 10~^ Mq year^^). 
Schoier et al 2002 assumed that the entire (~ 3000 AU in size) circumbinary envelope is 
infalling and derived the mass accretion rate from a model fitting to their JCMT spectra. 
On the other hand, from our higher spatial-resolution observations we suggest that only the 
compact HCN/HC^^N component associated with Source A is infalling, and we derived the 
mass accretion rate in the compact HCN component from our toy model. The difference of 
the estimated mass accretion rate probably arises from the different configuration adopted. 
From the estimated mass accretion rate, the accretion luminosity (= Lace) can be calculated 
as 



Lace — — 5 (3) 

where G is a gravitational constant, JVL is a mass of the central protostar, a nd i?* is 



a radius of the central protostar. If we assume i?* = 4 Rq (IStahler et al. 19801 ) . Lace is 
estimated to be ~ 1600 Lq with our M value and to be ~ 390 Lq with the M value by 
Schoier et al. (2002). These values are much higher than the bolometric luminosity of 116293 
(= 27 Lq; Walker et al. 1986), so that the so-called "luminosity problem" (Kenyon et al. 
1990) discussed by Ohashi et al. (1996) and Saito et al. (1996) for L1551 IRS5 also exists in 
116293. The reason of this discrepancy is still unclear. One possible explanation is that the 
infalling motion found in the present study is not onto the surface of the protostar itself, but 
onto the central rotationally-supported disk around the protostar. Here, the radius of the 
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terminal point in the gravitational potential is likely to be much larger than that adopted 
in the above estimates. If we assume the radius of the terminal point to be ~ 1 AU, our 
estimated mass infalling rate is consistent with the bolometric luminosity. Another possible 
explanation is that the accretion may be non-steady. 

We note that an infalling motion has the same radial dependence as that of a Keplerian 
rotation (oc r~^'^), and that they are indistinguishable from the examination of the P-V 
diagram alone. The only clue to distinguish the two motions is the disk orientation as 
compared to the direction of the velocity gradient. In fact, Huang et al. (2005) have 
interpreted the velocity gradient in the HC^^N emission as a Keplerian rotation around Source 
A. However, the position angle of the flattened structure in the optically- thin line (-16°) a s 



well as that of the inner rotating- infalling disk traced by the H2O maser (llmai et al. 19991 ). 
seems to favor the interpretation of the infalling motion. Furthermore, the interpretation 
of the Keplerian rotation implies that there is a well-developed 500 AU scale rotationally- 
supported disk around the Class protostar with the significant molecular outflows. In 
general, such well-developed large-scale (~ 500 AU) Keplerian disks are observed around 
more evolved sources, that is, Class II sourc es, as in the case of a circumbinary Keplerian 



ring arou nd GG Tau (IGuilloteau et al. 19991) and a Keplerian protoplanetary disk around 



DM Tau ( Guilloteau &: Dutrey 19981 ). From these considerations, we suggest that the infall 



gas motion is a more appropriate interpretation than that of the Keplerian rotation. 

The SMA+JCMT P-V diagram in Figure 7 exhibits another lower-velocity, extended 
component as seen in the velocity channel maps of Figure 4, which is mostly filtered out by 
the interferometer alone. The velocity structure in this extended lower- velocity component 
is clearly different from that in the compact high- velocity component; the velocity increases 
further from Source A. Since the cut of the P-V diagram is approximately parallel to the axis 
of the associated CO outflow (see Figure 6) and this type of velocity gradients is expected 
in molecular outflows driven by the wide-angle wind, we suggest that this extended HCN 
emission arises from turbulent gas swept up by the associated outflow from Source A. Our 
simple thin-layer model with the Gaussian intensity distribution shown in the right panel of 
Figure 7 is consistent with the interpretation that this low-velocity component could be an 
outflowing gas perpendicular to the compact high-velocity structure. 

In Figure 8, we show P-V diagrams along the cut through Source A and B (P. A. = 
138.8°) in the HCN emission with the SMA-only (gray) and SMA-^JCMT (black contour), 
and in the HC^^N emission with the SMA (red contour). At the position of Source A the wide 
velocity width in the HCN and HC^^N high- velocity component is evident. At the position of 
Source B, a weak HC^^N component with a narrow line width (~ 2 kms~^) is seen, which is 
smeared out in the total integrated intensity map of Figure 3 and embedded in the extended 
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HCN emission. This component at Source B is most likely a molecular structure found by 
the previous SMA observations (IKuan et al. 2004: IChandler et al. 20051) . as well as that by 
the previous Plateau de Bure observations (IBottinelli et al. 20041 ). In Figure 8, there is also 
an extended structure with narrow line width, which is likely to be the extended ambient 



gas. 



4.2. Different Evolution of the Binary Protostars ? 

Our submillimeter HCN and HC^^N observations of 116293 show that an intense compact 
(~ 500 AU) molecular component with a wide velocity width (> 10 kms~^) is associated 
with Source A but only a weak molecular component with a narrow velocity width (~ 2 
kms~^) is associated with Source B. SMA observations of 116293 in the CO (2-1, 3-2) 
lines have revealed that Source A is clea rly driving mole cular outflows, and that there is 



no clear outflow activity from Source B (lYeh et al. 20071 ). Source A is associated with a 



number of hot core trac ers such as complex organic molecules that are less prominent or 



not present at Source B (IKuan et al. 2004 : IChandler et al. 2005[ ) . These results imply that 



Source A and B are at different evolutionary stages in the common circumbinary envelope, 
although it is unclear which source is younger. In terms of protostellar activity traced by 
molecular lines. Source A seems to be younger than Source B, since Source A has substantial 
surrounding molecular gas and clear outf low emission, and shows a number of hot core tracers 



( iKuan et al. 2004 : IChandler et al. 2005 1). while Source B is more like a T Tauri star without 
a significant gas component or CO outflow. On the other hand, if the compactness in the 
continuum emission represents younger evolutionary stage of the circumstellar disk. Source 
B may be younger than Source A as suggested by Rodriguez et al. (2005) and Chandler 
et al (2005). In fact. Chandler et al (2005) have found redshifted SO (77-65) absorption 
against the strong continuum emission at Source B, which they suggest is the unambiguous 
detection of infall toward the central protostar. 

There are other examples of binary protostars where the components exhibit different 
observational characteristics that suggests they are in different evolutionary stages. Bourke 
(2001) has found two 7 sources with a separation of ~ 3400 AU in BHR 71 with the 
Infrared Space Observatory. The brighter companion (IRS 1) is associated with the millime- 
ter continuum source (BHR 71 mm) while the other weaker companion (IRS 2) is not, and 
IRS 1 drives a powerful molecular outflow while IRS 2 drives a much less massive outflow. 
From these results, it is suggested that IRSl and IRS2 in BHR 71 is a protobinary system 
at different evolutionary stages, and that IRS 2 is more evolved than IRS 1. High angular- 
resolution (~ 2") observations of CB 230 show compact millimeter emission associated with 
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only one component of an NIR protostellar pair separated by ~ 10", suggesting that like 
BHR 71 only this component has a substantial circumstellar disk, although both protostars 



in CB 230 drive CO outflows (iLaunhardt 200 ll ). A similar example is also seen in the SVS 



13 close binary system (separation ~ 65 AU). Anglada et al. (2004) reported from their 
high-resolution 7-mm observations with VLA that only one of the components of the SVS 
13 system (VLA 4B) is associated with detectable circumstellar dust emission, while the 
other component is optically visible without significant dust emission. These recent obser- 
vational results indicate that members of protostellar binary systems can exist at different 
evolutionary stages. 

Theoretically it is not clear how to form binary companions at different evolutionary 
stages in the common envelope and this idea remains controversial. According to recent the - 



oretical models of binary formation in this common envelope type ( INakamura fc Li 20031 ). 
fragmentation of "the first bar" which was made through the contraction of the envelope pro- 
duces binary or multiple stars, but the stellar age must be the same. Subsequent merging of 
the fragments may be necessary to make binary companions at different evolutionary stages. 
Ochi et al. (2005) have conducted high-resolution numerical simulations of an accretion from 
a circumbinary envelope onto the primary and secondary of the binary companion. Their 
simulations suggest that the accretion rate of the primary is larger than that of the secondary, 
regardless of the specific angular momentum of the accreting gas, and that the gas accretion 
tends to increase the mass difference between the primary and secondary. Their results are 
qualitatively different from earlier works by Bate & Bonnell (1997), which suggests that the 
primary accretes more than the secondary only when the accreting gas has a low specific 
angular momentum. More theoretical works, as well as an accumulation of observational 
data, are required to address this issue. 



5. Conclusions 

We have carried out submillimeter interferometric and single-dish observations of the 
well-known protobinary system IRAS 16293-2422 with the SMA in the HCN (4-3) and 
HC^^N (4-3) lines, in the continuum at 354.5 GHz and with the JCMT in the HCN (4-3) 
line. These data have provided us with the following results; 

1. The 354.5 GHz continuum observations with the SMA resolve the individual binary 
companions (Source A and Source B). The dust emission at Source A is elongated along the 
North-East to South- West direction (300 x 150 AU; P. A. = 33°), while that at Source B is 
more compact (150 x 140 AU). The elongation at Source A is likely to reflect the presence 
of the close binary system (Aa and Ab). The integrated submillimeter flux in both sources 
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is ~ 4.0 Jy. 

2. The optically-thin HC^^N (4-3) emission reveals a compact (~ 500 AU) flattened 
structure (P.A. = -16°) associated with Source A. There is a clear velocity gradient in this 
compact structure along the East (Blue) to West (Red) direction, which can be interpreted 
as an infalling gas motion in the flattened structure onto Source A with a central mass of ~ 
1 Mq and an inclination angle of ~ 57° from the plane of the sky. The combined HCN image 
shows an extended (~ 3000 AU) circumbinary envelope as well as the compact structure 
associated with Source A. In the extended circumbinary envelope there is also a low-velocity 
North-East (Blue) to South-West (Red) gradient, which may be due to gas swept up by the 
outflow. 

3. Contrary to Source A, at Source B there is only a weak molecular component with a 
much narrower line width (~ 2 kms~^), and no clear outflow activity. Thes e results, together 



with the fewer numb er of complex molecules associated with Source B (IKuan et al. 2004 



Chandler et al. 20051 ). may suggest different evolutionary stages between Source A and B in 
the common circumbinary envelope. With the data presented we are unable to determine 
whether Source B is younger or older. 
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A. Combining SMA and JCMT Data 

The procedure adopted to combine the SMA and JCMT data is similar to that described 
by Takakuwa et al. (2003a), which is based on the description for combining single-dish and 
interferometric data by Vogel et al. (1984), and on the MIRIAD scripts developed by Wilner 
& Welch (1994). We flrst re-sampled the JCMT image cube along the velocity axis to match 
the velocity channels with that of the SMA visibility. Next, we deconvolved the JCMT 
image by the 15" Gaussian image which represents the JCMT beam. Then, we multiplied 
the JCMT image with the 35" Gaussian image which is intended to approximate the SMA 
primary beam, that is, inverse of the primary beam correction. In these beam corrections 
we ignore any effect of sidelobes of the JCMT and SMA beam, which are not well-known. 
From the resultant JCMT image cube we made JCMT visibility data which flU in the central 
hole of the SMA uf-track by the MIRIAD task uvrandom and uvmodel. The single-dish 
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sampling in the uv plane was chosen so that the number of the uv data points per unit uv 
distance is similar to that of the SMA data (= 40 uv data points per IkX). We checked that 
the JCMT amplitude is approximately consistent with the SMA amplitude. 

The JCMT visibility data were Fourier-transformed along with the SMA visibility data 
to make a combined image cube by the MIRIAD task invert. We adopted robust = 0.5 
weighting for the imaging, which seems to provide a good compromise between the spatial 
resolution and the noise level. The resultant synthesized beam size and the rms noise level 
is 1'.'3 X 1'.'2 (P. A. = 30°) and ~ 0.90 Jy beam~^ channel"^, respectively. Negative sidelobes 
in the SMA synthesized beam are significantly suppressed in the combined beam. The final 
combined image with the JCMT and SMA has a high spatial dynamic range from large-scale 
(~ 3000 AU) to small scale (~ 200 AU) without any effect of missing flux. In Figure 9, 
we compare the JCMT, SMA, and combined images of 116293 in the HCN emission. The 
JCMT image shows a 3000 AU-scale circular blob and western elongation centered at Source 
A without any fine-scale structure. In the SMA image, finer structures are picked up but 
extended gas components are mostly resolved out. In the combined image, both extended 
and compact structures are successfully sampled. 



B. Importance of Combining Single-Dish and Interferometric Data 

Our combined SMA -|- JCMT observations of 116293 in the submillimeter HCN emission 
have revealed that the submillimeter emission extends more than ~ 3000 AU in 116293. This 
extent cannot be recovered with the SMA, since the minimum projected uv distance of the 
SMA observations was ~ 10.0 kX and the SMA observations were insensitive to structures 



more extended than ~ 2600 AU at 10% level (IWilner fc Welch 19941 ). In fact, our SMA 
observations recovered only ~ 35% of the total JCMT flux toward Source A. For comparison, 
SMA observations of L1551 IRS5, another low-mass protostellar envelope, recovered only ~ 
11% of the total submillimeter CS (7-6) flux, wh ich suggests that ther e is > 1500 AU-scale 



extended submillimeter emission in L1551 IRS5 (ITakakuwa et al. 20041 ) 



The rotational energy level of the upper state is 43 K and 66 K in HCN (4-3) and CS (7- 
6), and the critical density > 10^ cm~^ and > 10'' cm~^, respectively. These temperatures 
and gas densities are much higher than those traced by millimeter-wave tracers, such as 
C^^O (J=l-0) (5 K and ~ 10^ cm'^; Sargent et al. 1988, Momose et al. 1998) or H^3co+ 
(J=l-0) (4 K and ~ 10^ cm~^; Saito et al. 2001). It is puzzling why those submillimeter 
molecular lines which trace such warmer and denser gas can be so extended. It seems to be 
difficult to make gas temperature high enough only by the heating from central stars (e.g.. 
Lay et al. 1994). As discussed in §4.1, there is significant contamination from the associated 
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outflows in the submillimeter HCN emission, and the interaction between the circumbinary 
envelope and the outflowing gas could be a source of such extended submillimeter emission 
( lAvery &: Chiao 19961 : iNakamura 2000l ). To sample those extended submillimeter emission 
lines and to study the structure and kinematics from protostellar envelopes (> 1500 AU) to 
compact circumstellar disks (< 200 AU) unambiguously, it is quite important to make both 
submillimeter single-dish and interferometric observations and to combine the two data sets. 
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Fig. 1. — 354.5 GHz continuum image in 116293 taken with the SMA. Contour levels are 
from 3cr in steps of Sex (la = 0.066 Jy beam^^). The synthesized beam size, shown at the top 
left corner of the panel, is 1'.'05 x 0'.'64 (P.A. = 38.6°) with the uniform weighting. Crosses 
with the label "A" and "B" indicate positions of the 2.7-mm continuum peak observed with 
BIMA (Looney, Mundy, & Welch 2000), and those with the label "Aa" and "Ab" positions 
of the 300 GHz continuum peaks observed by Chandler et al. (2005). The size of the crosses 
shows the astrometric accuracy of the BIMA data. 
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Fig. 2. — HCN (4-3) line profiles toward Source A (left) and Source B (middle) taken with 
the SMA (bottom) and SMA + JCMT (middle), and HC^^N (4-3) line profiles taken with 
the SMA (top). For comparison with the JCMT HCN (4-3) spectrum, we also plot the 
average spectra toward Source A and Source B (right). The spatial resolution of the HCN 
spectra is 1'.'3 x 1".2 (P. A. = 30°). The spatial resolution of the JCMT spectra is 15", and 
that of the SMA HC^^N spectra 1'.'6 x 1'.'3 (P.A. = 14°). The rms noise level of the JCMT, 
SMA, and the combined spectra is 0.13 K, 7.4 K, and 5.6 K, respectively. 
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Fig. 3. — Total integrated intensity maps in 116293 taken with the SMA (left) and 
SMA+JCMT (middle) in the HCN (4-3) emission, and the HC^^N (4-3) emission with the 
SMA (right). Crosses in each panel indicate positions of the 354.5 GHz continuum emission 
observed with the SMA (Fig.l, Table 2). Filled ellipses at the top left corner of each panel 
show the synthesized beam, that is, 1'.'3 x 1'.'2 (P.A. = 30°) in the SMA and the combined 
map, and 1'.'6 x 1'.'3 (P.A. = 14°) in the SMA HC^^N map. Contour levels are from 33.50 
K km s~^ in steps of 22.33 K km s~^ in the SMA and combined HCN map {3a in steps of 
2a in the SMA HCN map), and from 3a in steps of 2a {la = 6.49 K km s'^) in the SMA 
HC^^N map. 
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Fig. 4.— SMA + JCMT velocity channel maps of 116293 in the HCN (4-3) hne. Each 
panel is an average of 3 channels (1 channel = 0.172 km s~^). Crosses in each panel indicate 
positions of the 354.5 GHz continuum emission observed with the SMA (Fig.l, Table 2). 
The synthesized beam size is shown as a filled ellipse at the top left corner of each panel. 
Numbers at the upper-left corner of each panel are Vlsr (km s~^). Contour levels are from 
2(7 in steps of 2a {la ^ 3.23 K). 
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Fig. 5.— SMA + JCMT velocity channel maps of 116293 in the HCN (4-3) emission at 
four different velocity ranges: (a) low-velocity blueshifted range (2.54 - 3.40 km s~^); (b) 
low-velocity redshifted range (4.95 - 6.50 km s^^); (c) high-velocity blueshifted range (-2.62 
- -1-0.65 km s~^); and (d) high-velocity redshifted range (7.36 - 9.08 km s~^), superposed on 
the CO (2-1) outflow map taken with the SMA (blue and red contours; Yeh et al. 2007). 
Contour levels are from 5.30 (K) in steps of 3.54 (K) in the HCN map and 5.68 (K) in steps 
of 5.68 (K) in the CO map. Green crosses indicate positions of the protobinary (Fig.l, Table 
2). A filled ellipse at the top left corner of panel (a) indicates the synthesized beam of the 
HCN observations (1'.'3 x 1'.'2; P.A. = 30°), while an open ellipse at the top left corner of 
panel (b) the synthesized beam size of the CO map (3'.'3 x 2'.'0; P.A. = 40°). 
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Fig. 6. — High-velocity blueshifted {Vlsr = -2.62 - 2.37 km s ^) and redshifted (Vlsr = 5.98 
- 9.08 km s-i) HCN (4-3) emission taken with the SMA and JCMT (left) and HC^^N (4-3) 
emission taken with the SMA (right) in 116293. The velocity range is set to be wider than 
that of Figure 5 in order to improve the signal-to-noise ratio of the HC^^N image. The CO 
(2-1) outflow map taken with the SMA (Yeh et al. 2007) is overlaid on the HCN and HC^^N 
images. The blueshifted and redshifted components in the HCN and HC^^N lines are shown 
in blue and red contours, while the blueshifted and redshifted CO outflow in black bold and 
thin contours, respectively. Contour levels are from 19.4 (K) in steps of 3.5 (K) in the HCN 
map, from 3.2 (K) in steps of 2.1 (K) in the HC^^N map, and from 17.1 K in steps of 11.4 K 
in the CO map. Crosses indicate positions of the 354.5 GHz continuum emission observed 
with the SMA (Fig.l, Table 2). A filled ellipse at the top left corner shows the synthesized 
beam of the HCN (1'.'3 x 1'.'2; P.A. = 30°) and HC^^N (1'.'6 x 1'.'3; P.A. = 14°) images, while 
open ellipses at the bottom left corner the synthesized beam of the CO (2-1) image (3'.'3 x 
2'.'0; P.A. = 40°). Lines a-b and c-d indicate cuts of the Position- Velocity diagrams shown 
in Figure 7 and 8, respectively. 
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Fig. 7. — Left: Position- Velocity (P-V) diagrams along the minor axis of the compact 
flattened structure at Source A in the HC^^N (4-3) emission (red contour; P. A. = 74°), as 
well as in the HCN (4-3) line with the SMA + JCMT (black contour) and SMA (gray). 
Contour levels are from 11.2 (K) in steps of 11.2 (K) in the HCN emission and from 4.23 (K) 
in steps of 4.23 (K) in the HC^^N emission. A horizontal dashed line indicates the position 
of Source A, and a vertical dashed line the systemic velocity (~ 3.6 km s~^). There are two 
distinct components, one is the high-velocity compact structure seen in the SMA HCN and 
HC^^N images and the other the extended low-velocity envelope seen in the combined image. 
Right: Model P-V diagrams to interpret the observational P-V diagrams in the left panel. 
Contour levels are the same as in the combined HCN P-V diagram. The P-V diagram of the 
high-velocity compact structure can be interpreted as an accretion motion onto the central 
mass of 1 M0 in the geometrically thin disk with an inclination angle of 57° from the plan 
of the sky and an internal velocity dispersion (= a) of 1.0 km s~^. The P-V diagram of the 
extended component can be understood as an outflowing gas perpendicular to the flattened 
structure with an internal velocity dispersion of 0.3 km s~^. The model emission is assumed 
to have Gaussian intensity distribution. The strong self-absorption around Vlsr = 4-5 km 
s~^ seen in the real data can not be reproduced with this "toy" model. 
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Fig. 8.— P-V diagrams in the HCN (4-3) line with the SMA + JCMT (black contour) and 
SMA (gray) and in the HC^^N (4-3) line with the SMA (red contour), along the axis joining 
sources A and B (P.A. = 139°). Contour levels are from 11.2 (K) in steps of 11.2 (K) in the 
HCN emission and from 4.23 (K) in steps of 4.23 (K) in the HC^^N emission. Upper and 
lower horizontal dashed lines show the position of Source A and Source B, respectively, and 
a vertical dashed line the systemic velocity (~ 3.6 km s~^). A solid horizontal line indicates 
the middle point between Source A and Source B, and should be the binary axis on the 
assumption of the same mass in Source A and Source B. 
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Fig. 9. — Total integrated intensity maps in 116293 taken with JCMT (upper), the SMA 
(lower-left), and SMA+JCMT (lower-right) in the HCN (4-3) emission. An open circle in 
the JCMT map indicates the field of view of the SMA, and crosses in each panel indicate 
positions of the protobinary (Fig.l, Table 2). Filled ellipses at the top left corner of each 
panel show the synthesized beam, that is, 15" in the JCMT map and 1'.'3 x 1'.'2 (P.A. = 30°) 
in the SMA and the combined map. Contour levels are from 4.76 K km s~^ in steps of 4.76 
K km s~^ in the JCMT map, and from 33.50 K km s~^ in steps of 22.33 K km s~^ in the 
SMA and combined map. 



-30- 



Table 1. Parameters for the SMA Observations of IRAS 16293-2422 



Parameter 




Value 


2003 March 14 


2003 July 12 2004 June 19 


Number of Antennas 


5 


4 5 


Right Ascension (J2000) 




le'' 32™ 22^91 


Declination (J2000) 




-24°28' 35^52 


Primary Beam HPBW 




~ 35" 


Synthesized Beam HPBW (HCN) 


l'.'3 


X 1^'2 (P. A. = 30°) 


Synthesized Beam HPBW (HCi'^N) 


1"6 


X 1{'3 (P. A. = 14°) 


Synthesized Beam HPBW (Continuum; uniform) 


I'.'l 


X O'.'e (P.A. = 39°) 


Baseline Coverage 




10 - 223 (kA) 


Conversion Factor (Line) 


1 (Jy beam-i) = 6.2 (K) 


Frequency Resolution (HCN) 


203.125 kHz ~ 0.172 km s'^ 


Frequency Resolution (HC^^N) 


812.5 kHz ~ 0.708 km s'^ 


Bandwidth 


82 MHz X 8 


82 MHz X 12 82 MHz x 24 


Gain Calibrator 


nrao530 


1743-038 1743-038 


Flux of the Gain Calibrator 


1.5 Jy 


1.5 Jy 2.3 ,ly 


Passband Calibrator 


Mars 


Mars, Uranus Juijiter, Uranus 






3c279, 1924-292 


Flux Calibrator 


Callisto 


Uranus Uranus 


System Temperature (DSB) 




~ 350 K 


rms noise level (Continuum) 




' 0.06 Jy beam-i 


rms noise level (Line; SMA + JCMT) 


~ 0.9 Jy beam-i / 203.125 kHz 
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Table 2. 354.5 GHz Continuum Results in 116293 with the SMA 



Protostar 


R.A.^ 




Decl.'' 


Deconvolved FWHM Size"" 


P.A.=- 


Total Flux=- 




(J2000) 




{J2000) 


(AU X AU) 


(°) 


(mjy) 


Source A 




.87 


-24°28'36'.'6 


300 X 150 


32.6 


3840 


Source B 


16>'32™22= 


.62 


-24°28'32'.'6 


150 X 140 


24.8 


4050 



''Estimated from 2-dimensional Gaussian fittings to the image. 



